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Transport coefficients associated with the mass flux of a binary mixture of Maxwell molecules 
under uniform shear flow are exactly determined from the Boltzmann kinetic equation. A normal 
solution is obtained via a Chapman-Enskog-like expansion around a local shear flow distribution 
- that retains all the hydrodynamics orders in the shear rate. In the first order of the expansion 

I the mass flux is proportional to the gradients of mole fraction, pressure, and temperature but, 

^— ■») ■ due to the anisotropy induced in the system by the shear flow, mutual diffusion, pressure diffusion 
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and thermal diffusion tensors are identified instead of the conventional scalar coefficients. These 
tensors are obtained in terms of the shear rate and the parameters of the mixture (particle masses, 
concentrations, and force constants). The description is made both in the absence and in the 
presence of an external thermostat introduced in computer simulations to compensate for the viscous 
heating. As expected, the analysis shows that there is not a simple relationship between the results 
, with and without the thermostat. The dependence of the three diffusion tensors on the shear rate is 

illustrated in the tracer limit case, the results showing that the deviation of the generalized transport 
coefficients from their equilibrium forms is in general quite important. Finally, the generalized 
transport coefficients associated with the momentum and heat transport are evaluated from a model 
kinetic equation of the Boltzmann equation. 
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PACS numbers: 51.10.+y, 05.20.Dd, 05.60.-k, 47.50.-d 



! I- INTRODUCTION 

■ 

^ , The description of transport properties for states close to equilibrium in gaseous binary mixtures is well established. 
Ch In these situations, the Curie principle^ states that the presence of a velocity gradient (second-rank tensorial quantity) 
, cannot modify a vectorial quantity such as the mass flux ji, which is generated by gradients of mole fraction xi, 
(-H pressure p, and temperature T. As a consequence, the mutual diffusion coefficient D (which couples the mass current 
Q I with Vxi), the pressure diffusion coefficient Dp (which couples the mass current with Vp) and the thermal diffusion 
O 1 coefficient Dt (which couples the mass current with VT) do not depend on the velocity gradient. However, when the 
shear rate applied is large, non-Newtonian eff'ects are important so that the Curie principle does not hold and the 
, coefficients associated with the mass transport are affected by the presence of shear flow. In particular, if the spatial 
^ ' gradients Vxi, Vp, and VT are weak, one expects that the flux ji is still linear in these gradients but the standard 
' scalar coefficients {D, Dp, Dt} must be replaced by the shear-rate dependent second-rank tensors {Dij,Dp^ij,DT,ij}- 
' The aim of this paper is to determine the above tensors in the framework of the Boltzmann equation. 
I We are interested in a situation where weak spatial gradients of mole fraction, pressure, and temperature coexist 
. . with a strong shear rate. Under these conditions, the application of the conventional Chapman-Enskog expansion^ 
On around the local equilibrium state to get higher order hydrodynamic effects (Burnett, super-Burnett, . . .) to the mass 
flux turns out to be extremely difficult. This gives rise to look for alternative approaches. A possibility is to expand 
around a more relevant reference state than local equilibrium. Since we want to compute the mass transport in a 
strongly sheared mixture, the so-called uniform shear flow (USF) state can be chosen as the reference state. The USF 
state is characterized by constant mole fractions, a uniform temperature, and a linear velocity profile Ux — ay, where 
a is the constant shear rate. Due to its simplicity, this state has been widely used in the past to shed light on the 
complexities associated with the nonlinear response of the system to the action of strong shearing. In addition, the 
' USF state is one of the rare exceptions for which the hierarchy of moments of the Boltzmann equation admits an 
exact solution for single^ and multicomponent gases^ of Maxwell molecules (repulsive potential of the form r~^). In 
this case, explicit expressions of the pressure tensor (which is the relevant irreversible flux of the problem) have been 
obtained for arbitrary values of the shear rate and the parameters of the system (masses, concentrations and force 
constants) . 

As said before, here we want to compute the mass transport under USF for Maxwell molecules. Since the mixture 
is slightly perturbed from the USF, the Boltzmann equation can be solved by an expansion in small gradients around 
the (local) shear flow distribution instead of the (local) equilibrium. This is the main feature of the expansion since 
the reference state is not restricted to small values of the shear rate. In the flrst order of the expansion, the set 
of generalized transport coefficients {Dij, Dp^ij, DT,ij} are identified from the mass flux ji as nonlinear functions of 
the shear rate and the parameters of the mixture. This Chapman-Enskog-like expansion has been used to analyze 
transport properties in spatially inhomogeneous states near USF in the case of ordinary gases^i^ and more recently 
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in the context of granular gasesT^^S, 

Some previous attempts have been carried out earher by the author and coworkersi^iiiiiSiiHiii in the case of the 
diffusion tensor Dij. However, aU these studies have been restricted to perturbed steady states with the constraints 
p = const and T — const. Ahhough steady states are in general desirable for practical purposes, especially in computer 
simulations,^^ here we extend the above studies to a general time and space dependence of the hydrodynamic fields. 
This allows us to evaluate new contributions to the mass flux (those proportional to Vp and VT), which where not 
taken into account in the previous studies;^ 

The plan of the paper is as follows. First, a brief summary of the results obtained from the Boltzmann equation 
for a binary mixture of Maxwell molecules under USF is presented in Sec.|TTl Section Hill deals with the perturbation 
scheme used to solve the Boltzmann equation for the mixture to first order in the deviations of the hydrodynamic 
field gradients from their values in the reference shear flow state. The generalized transport coefficients characterizing 
the mass transport around USF are also defined in Sec. IIIII These coefficients are explicitly obtained in Sec. IIVI with 
and without the presence of an external thermostat introduced usually in computer simulations to compensate for 
the viscous heating. The dependence of some of these coefficients on the shear rate is illustrated with detail in the 
tracer limit case, showing that the influence of shear flow on mass transport is quite significant. The paper is closed 
by a brief discussion of the results in Sec. |Vl the generalized transport coefficients associated with the momentum 
and heat transport are evaluated in Appendix [D] from a simple model kinetic equation of the Boltzmann equation. 



II. A BINARY MIXTURE UNDER UNIFORM SHEAR FLOW 



We consider a dilute binary mixture where /s(r,v;t) is the one-particle velocity distribution function of species s 
(s = 1, 2). The time evolution of the distributions fg is given by the set of two coupled nonlinear Boltzmann equations: 



dt+^-^ + j. ^ /,(r,v,i)=^J,,[v|/4i),/r(i)] , (2.1) 

^ r— 1 

where rris is the mass of a particle of species s, F^ is a possible external force acting on particles of species s, and 
Jsr \y\fsi fr] is the Boltzmann collision operator, which in standard notation reads^ 



JsAfsJr] ^ j '^''^j C«^^|V- Vl|(7,,(v- Vi,0)[/,(v')/.(v;) -/,(v)/,(vi)] . (2.2) 

The basic moments of fa are the species number densities 

j dv/„ (2.3) 

and the mean velocity of species s 

Us^ — f dvv/,. (2.4) 

ris J 

These quantities define the total number density n = ^^"^ ^he flow velocity u = Ps'^s/p, where ps = msUs 

is the mass density of species s and p = Ps the total mass density. The temperature T is deflned as 



where /cb is the Boltzmann constant and V = v — u is the peculiar velocity. The second identity in ()2.5p defines the 
partial kinetic temperatures of species s. They measure the mean kinetic energy of particles of species s. Moreover, 
in a dilute gas the hydrostatic pressure p is given hy p = nkBT. The quantities Us, u, and T are associated with 
the densities of conserved quantities (mass of each species, total momentum, and total energy). The corresponding 
balance equations define the dissipative fluxes of mass 



=m,ydvV/„ (2.6) 

momentum (pressure tensor) , 

P = 5] Ps=5] m.ydvVV/,, (2.7) 
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and energy (heat flux) 



(2.8) 



The second equahties in Eqs. ()2.7p and (j2.8|) define the partial contributions and to the pressure tensor and 
heat flux, respectively. The fact that the mass flux is defined with respect to the local center-of-mass velocity u 
implies that 



E 



0. 



(2.9) 



The USF state is macroscopically defined by constant densities rig, a spatially uniform temperature T(t) and a 
linear velocity profile u{y) = Ui(y) = U2(y) = ayx, where a is the constant shear rate. Since Us and T are uniform, 
then = q = 0, and the transport of momentum (measured by the pressure tensor) is the relevant phenomenon. 
In the USF problem, the temperature tends to increase in time due to viscous heating. Usually, an external force 
(thermostat) is introduced in computer simulations to remove this heating effect and keep the temperature constantf^^ 
The simplest choice is a Gaussian isokinetic thermostat given by 



(2.10) 



where the thermostat parameter a is a function of the shear rate adjusted as to keep the temperature constant. 
The implicit assumption behind the introduction of these forces is that they play a neutral role in the transport 
properties, so that the latter are the same with and without a thermostat, when conveniently scaled with the thermal 
speed. Nevertheless, this expectation is not in general true, except for some specific situations and/or interaction 
potentials i^ii^ 

At a microscopic level, the USF is characterized by a velocity distribution function that becomes uniform in the 
local Lagrangian frame, i.e., /s(r,v;t) = /s(V,t). In that case, Eq. (|2.ip with the choice (|2.10|) reduces to 



Jll[flJl]+Mflj2] 



(2.11) 



and a similar equation for /2. Here, ay — aSixSjy. The hierarchy of velocity moments associated with the Boltzmann 
equation (|2.11[) can be recursively solved in the particular case of Maxwell molecules, i.e., when particles of species r 
and s interact through a potential of the form Vrs{r) = Krsr~'^. The key point is that for this interaction the collision 
rate garsig, 0) is independent of the relative velocity g and so the coUisional moments of order k only involve moments 
of degree smaller than or equal to k. In particular, the first- and second-degree coUisional moments are given b y^^i^^ 



m. 



J dvY, Jsrif s, fr 



A. 



msruir 



{Ps'ir - PtU) , 



(2.12) 



rUs J dvWJsrifsJr] 



{nia + mr)m, 

-psPr - PrPs + jjr + jrjs] 



PsPr + PrPs - gjs • jr ) 1 



{nis + TOr)ms 



2( — p,P, -p,P, 



(jjr + jrjs) 



where Ps — -^trP^ = UgkBTs is the partial hydrostatic pressure and 



A,r — 1.697r K 



iris + ITT-r 



1/2 



A' = 2.6l7r K 



rrig + rrir 



1/2 



(2.13) 



(2.14) 



Thanks to the above property, exact expressions of the pressure tensor P for a binary mixture of Maxwell molecules 
under USF were obtained some time agO;^ The nonzero elements of P are related to the rheological properties of the 
mixture, namely, the nonlinear shear viscosity and the viscometric functions. In reduced units, they turn out to be 
nonlinear functions of the (reduced) shear rate a* = a/C, (where C is a convenient time unit defined below) and the 
parameters of the mixture: the mass ratio p = mi/m2, the mole fraction xi — ni/n and the force constant ratios 
^11/^12 and K22/ 1^12- It must be noted that in the particular case of Maxwell molecules there is an exact equivalence 
between the USF results with and without the external forces (|2.10p . As will be shown below, beyond the USF 
problem, the presence of the thermostat does not play a neutral role in the results and a certain infiuence may exist. 
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III. CHAPMAN-ENSKOG-LIKE EXPANSION AROUND USE 



As said in the Introduction, the main aim of this work is to analyze mass transport of a dihite binary mixture 
subjected to USF. In that case, let us assume that the USF state is disturbed by small spatial perturbations. The 
response of the system to those perturbations gives rise to contributions to the mass flux that can be characterized 
by generalized transport coefficients. This Section is devoted to the evaluation of those coefficients. 

In order to analyze this problem we have to start from the set of Boltzmann equations (|2.ip with a general time and 
space dependence. Let Ug = a • r be the flow velocity of the undisturbed USF state, where the elements of the tensor a 
are = aSixSjy. In the disturbed state, however the true velocity u is in general different from Uq, i.e., u = Uq + Su, 
du being a small perturbation to Uq. As a consequence, the true peculiar velocity is now c = v — u = V — 6u, where 
V = V — Uq. In the Lagrangian frame moving with ug, the Boltzmann equations (12. 1[) can be written as 



^/i - ^ {a^JV, + aV,) /i + (V + uo) ■ V/i + aSu ■ ^ = Jii[/i, /i] + Ji2[/i, M, (3.1a) 



- ^ (oyl/, + aV,) /2 + (V + Uo) • V/2 + adu-^ = Jsai/a, M + J2i[/2, /i], (3.1b) 



where here the derivative V/s is taken at constant V. In addition, in Eqs. (13. lap and (I3.1bp the thermostat force has 
been assumed to be proportional to the actual peculiar velocity, Fg = —msa{Y — 6u) where now the parameter a 
is in general a function of r and t through their functional dependence on the hydrodynamic fields Us and T. The 
generalization of a to the inhomogeneous case is essentially a matter of choice. Here, for the sake of simplicity, we 
will take two different choices for a: (i) a = 0, so that the temperature grows in time, and (ii) the same expression 
obtained in the (pure) USF problem, except that the densities and temperature are replaced by those of the general 
inhomogeneous state. 

The macroscopic balance equations associated with this disturbed USF state are obtained by taking moments in 
Eqs. ([3Ta|) and (IS.lbj) with the result 

dtn, + Uo • Vn, + V • (nju) = (3.2) 

rris 

dtSui + aijSuj + (uo + Su) ■ VSui = —p^^\/jPij, (3.3) 

2 

^ndtT + |n(uo + 5u) ■ VT = -aP,y + |t ^ - ( V • q + P : VJu + 3pa) , (3.4) 

where the mass flux j^, the pressure tensor P, and the heat flux q are defined by Eqs. (|2.6p . (|2.7p . and (|2.8p . respectively, 
with the replacement V ^ c. The corresponding balance equations for the mole fraction x\ — ni/n and the pressure 
p — nksT can be obtained from Eqs. (|3.2p and (|3.4p . They are given by 

(9txi + (uo + (5u) • Vxi = V-ji, (3.5) 



dtp + (uo + (5u) • Vp + pV • (5u = -- {aP^y + V • q + P : V^u + 3pa) . (3.6) 



We assume that the deviations from the USF state are small. This means that the spatial gradients of the hy- 
drodynamic fields are small. For systems near equilibrium, the specific set of gradients contributing to each flux is 
restricted by fluid symmetry, Onsager relations, and the form of entropy production.^ However, in far from equilibrium 
situations (such as the one considered in this paper), only fluid symmetry applies and so there is more flexibility in 
the representation of the heat and mass fluxes since they can be defined in a variety of equivalent ways depending 
on the choice of hydrodynamic gradients used. In fact, some care is required in comparing transport coefficients in 
different representations using different independent gradients for the driving forces. Here, as in previous works 
the mole fraction xi, the pressure p, the temperature T, and the local flow velocity Su are chosen as hydrodynamic 
fields. 

Since the system is strongly sheared, a solution to the set of Boltzmann equations (|3.1ap and (jS.lbp can be obtained 
by means of a generalization of the conventional Chapman-Enskog method^ in which the velocity distribution function 
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is expanded around a local shear flow reference state in terms of the smah spatial gradients of the hydrodynamic 
fields relative to those of USF. This is the main new ingredient of the expansion. This type of Chapman-Enskog- 
like expansion has been already considered to get the set of shear-rate dependent transport coefficienta^^ii^ in 
thermostatted shear flow problems and it has also been recently used for inelastic gaseSf^i^ii^ 
In the context of the Chapman-Enskog method,^ we look for a normal solution of the form 

/,(r,V,i) = /,[A(r,t),V], (3.7) 

where 

A(r,i) = {a;i(r,i),p(r,t),r(r,t),Mr,t)}- (3.8) 

This special solution expresses the fact that the space dependence of the reference shear flow is completely absorbed 
in the relative velocity V and all other space and time dependence occurs entirely through a functional dependence 
on the fields A{r, t). The functional dependence (|3.4p can be made local by an expansion of the distribution function 
in powers of the hydrodynamic gradients: 

/. [A{r, t , V] = /i") (V) + (V) + . . . , (3.9) 

where the reference zeroth-order distribution function corresponds to the USF distribution function but taking into 
account the local dependence of the concentration, pressure and temperature and the change V ^ V — 5u{r, t) = c. 

The successive approximations fi are of order k in the gradients of xi, p, T, and ou but retain all the orders in the 
shear rate a. Here, only the first-order approximation will be analyzed. 

When the expansion (|3.9p is substituted into the definitions (|2.6p . (|2.7p . and (|2.8p . one gets the corresponding 
expansions for the fluxes: 

j,s=ji°^+ji') + --- , (3.10a) 

p = p(") + p(i) + . . . , q = q(") + q(i) + . . . . (3.10b) 

Finally, as in the usual Chapman-Enskog method, the time derivative is also expanded as 

9t = c>i°'+9i'^+af^ + --- , (3.11) 

where the action of each operator cJ^'^'' is obtained from the hydrodynamic equations (j3.2p - (|3.4p . These results provide 
the basis for generating the Chapman-Enskog solution to the Boltzmann equations p. lap and p.lbp . 



Zeroth-order approximation 



Substituting the expansions (|3.10ap - (|3.1ip into Eq. (|3.1ap . the kinetic equation for /{ ' is given by 

|/^°^ - ^ ia.,V, + aV.) /}°) + (V + uo) • V/'"' + -Su ■ ^ = Jn[/f \ /f ^] + MA'\ A% (3-12) 
To lowest order in the expansion the conservation laws give 

afxi^O, T-ia^")T = p-ia(%=-A«p(0)_2a, (3.13) 

dl°^6u, + a,jSuj ^0. (3.14) 

If a = 0, then T-^dl"^T = p-^^f V = -2aP^°y^ /3p while if a = -aPi,°y^ /3p then 9^ ¥ = ^ = 0. 
Since /[^"^ is a normal solution, the time derivative in Eq. (|3.12p can be represented more usefully as 

T-f(O) T-f(O) »-f(0) af(0) 

^1 " dx, ""'^ dp dT ^ ^ dSui' ^ ' 



(3.15) 
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where in the last step we have taken into account that /^"^ depends on Su only through the peculiar velocity c. 
Substituting Eq. (j3.15|) into Eq. ()3.12|) yields the following kinetic equation for /}°'': 



(0) Mh 



(3.16) 



A similar equation holds for f2^\ The partial pressure tensors P^'^'' and can be obtained from Eq. (|3.16p and its 

counterpart for /2'''' when one multiplies both equations by m^cc and integrate over c. Their explicit forms can be 
found in the Appendix of Ref. 0. 



(0) d(o) 



B. First-order approximation 

The analysis to first order in the gradients is worked out in Appendix lAl The distribution function /|^' is of the 
form 



/(I) = Ai ■ Vxi + Bi • Vp + Ci • VT + Pi : V5u, 



(3.17) 



where the vectors {^1, J3i, Ci}, and the tensor Vi are functions of the true peculiar velocity c. They are the solutions 
of the following set of linear integral equations: 



— aPi°) + 2a j {pdp + TOt) Ai 



d d 



'^^y-g^ + -c] Ai+ CiAi + M1A2 



(3.18) 



—aPi°^ + 2a ) {pdp + Tdr) Bi 



acy-g^ + a— -0)61+ CiBi + M1B2 



2a 



-9pPi°)+ 2(1 



'^-^[l-pdp)Pi^) -2Tdpa 



(3.19) 



^aPi°J + 2o^ {pdp + Tdr) Ci - [acy^ + « A . c ) Ci + C^C^ + M1C2 



2a 
3p 



(1 + TdT)P^'(l + 2(1 + TdT)a 



2a 



dTP^f+2pdTa) Bi, 



(3.20) 



Pi°J+2a) {pdp + TdT)Vi,i 



3p 



d d 



(3.21) 



where Ai(c), Bi(c), Ci(c), and Di(c) are defined by Eqs. (|A8p - (|All[) . respectively. In addition, Ci and Aii are the 
linearized Boltzmann collision operators around the reference USE state: 



CiX = - [Jii[f['\x] + Jii[X,/f'] + Ji2[X,/f ]) 



(3.22a) 



(3.22b) 



In this paper we are mainly interested in evaluating the first-order contribution to the mass fiux It is defined 



as 



(3.23) 
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Use of Eq. (|3.17p into Eq. (|3.23p gives the expression 



.(1) mirn^n dxi p dp p dT 



where 

D 



27 — 

777.271 



y dcci A,i(c), (3.25) 



rriip 
P 



(3.26) 



^ y dcci Cij (c). 



miT 
P 



(3.27) 



Upon writing Eqs. p.25p - p.27p use has been made of the symmetry properties of ,4-1, Bi, and Ci. In general, the set 
of generalized transport coefficients Dij , , and Dt^zj are nonUnear functions of the shear rate and the parameters 
of the mixture. It is apparent that the anisotropy induced by the presence of shear flow gives rise to new transport 
coefficients for the mass flux, reflecting broken symmetry. According to Eq. (|3.24p , the mass flux is expressed in terms 
of a diffusion tensor Dij, a pressure diffusion tensor DpAj, and a thermal diffusion tensor DT,ij- 

To get the explicit dependence of the above transport coefficients on the parameter space of the problem, the form 
of a must be chosen. As said before, two choices will be considered here: (i) the unthermostatted case a — Q, and (ii) 

the thermostatted case a = —aPj^^y/ip. Both cases will be separately studied in the next Section. 

IV. MASS TRANSPORT UNDER SHEAR FLOW 

This Section is devoted to the determination of the generalized transport coefficients D^ , -D^.y , and D^^ij associated 
with the mass transport for the two choices of the thermostat parameter. These coefficients are given in terms of the 
solutions to the integral equations (|3.18p ~ (|3.20p . 



Unthermostatted USF state 



In the absence of an external thermostat (a = 0), the integral equations ()3.18p - ()3.20p become 

op OCx 

= Ai + g(pei + rCi)(9,,pi;;)), 



(4.1) 



-^aPi;;) {pdp + TOt) Bi - (^y9pPi°) + ac,^ ) Bi + CiBi + M1B2 



mj 

2aT 



Bi-^Ci(l-p9,)Pi°), 



(4.2) 



-^aPi°) {pdp + Tdr) Ci - ^(1 + raT)Pi°) + acy ^ 
= Ci + ^Bi(9TPi°)). 



9ct 



Ci + LiCi + Mie2 



(4.3) 



The dependence of P^j"^ on the pressure p and temperature T occurs explicitly and through its dependence on the 
reduced shear rate a* = a/C^. Here, the effective collision frequency C is given by^ 



C = 277 



p A'i2 



7771 + 7772 ksT 7771 + 7772 



(4.4) 
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where X'12 is defined in Eq. (|2.14p . Consequently, 



dpP^f = dppP:,ia*) = ( 1 - a*^ 1 (a*), (4.5) 



d 



(4.6) 



where P*j = P^j^ /p. In addition, the dependence of P^p on the mole fraction xi is also rather intricate and so the 
derivatives with respect to xi must be carried out with care.^^ The generalized coefficients Dij, Dp^tj, and DT,ij can 
be obtained from Eqs. (|4.ip - (|4.3p when one multiplies those equations by miCi and integrates over c. After some 
algebra (some technical details are provided in Appendix |B]) , one arrives at the following set of coupled algebraic 
equations: 



pX 



12 



TOim2 



-n P* 



5ik + ai, 



mim2 



B P* 



Pi 



8 P* 



3mim2n 



(4.7) 



- ^(1 - a*a,.)^;j '^.fc + a.fc 
mim2 op " ' 



DpM, - ^^{l-a*da.)^Pl,,-^P* 



2a 



a*DT,^Ji^a^P:y), 



(4.8) 



2a 



\m1m2 



{l + a*da')P*]S,k + 



O-ik 



p 



Pi 



DT,kj = -a*da' I P^ 

+ '^a*D,,,{da^P:^) 



p* 



(4.9) 



where P* — P'^s^ /p and use has been made of the relations (|4.5p and (|4.6p . 

In the absence of shear field (a = 0), then P^ij = XgSij, and P*j = 6ij, so that Eqs. (I4.7p - ()4.9p have the solutions 
Dij = Df)Sij, Dp ij = Dp oSij, and -Dr.i 



coefficients for Maxwell molecules <^ Their expressions are 



Ai 



0, where Dq and Dp Q are the conventional Navier-Stokes transport 

(4.10) 



Dp,a = — 3-("^2 - mi)Do. 



In this case the mass flux j'^' can be written asi^ 



Ji 



where 



mim2PiP2 „ (V<^i)t - (V</>2)i 

— ; — 2 — ^0 ^ 

ksp^ T 
T 



= — Vln(xsp), 



(4.11) 



(4.12) 



(j)s being the chemical potential per unit mass. The fact that the thermal diffusion coefficient vanishes when a* — Q 
is due to the interaction potential considered (Maxwell molecules) since this coefficient is different from zero for more 
general interaction potentialsi^ However, when the mixture is strongly sheared, the Boltzmann equation leads to 
contributions to the mass flux proportional to the thermal gradient, even for Maxwell molecules. 

p(0) 



and so Dr, 



0. Moreover, Eq. (|4.7p reduces to 
-1 



In the case of mechanically equivalent particles (/i — 1, kh = K22 — ^12), P*^ij = xiP*j, dxiPfJj — Pi'ij/xi 



AO) 



D 



O-ik 



xy \ ' n\i2lm - laP*y 



p: 



(0) 



kj 



(4.13) 



Equation (|4.13p is consistent with previous results derived for the self-diffusion tensor Furthermore, known results 
for the diffusion tensor^*' are also recovered in the tracer limit (xi —> 0). 



9 




10 



0.0 


-=:^= u=0.2 - 


-0.5 








-1.0 








-1.5 




M=JO\ 




-2.0 




1,1, 





0.0 0.4 0.8 1.2 1.6 2.0 



a* 

FIG. 3: Shear-rate dependence of the trace l-Dr.fcfc oi the thermal diffusion tensor for xi = 0, K22 = K12 and several values of 
the mass ratio fi = m\/m2. 




a* 



FIG. 4: Shear-rate dependence of the off-diagonal element —D*y of the mutual diffusion tensor for xi = 0, K22 = /ti2 and 
several values of the mass ratio /^ = mi/m2. 




a* 



FIG. 5: Slicar-ratc dependence of the off-diagonal element —D^^^y of the pressure diffusion tensor for xi = 0, K22 = K12 and 
several values of the mass ratio /x = mi/ 1712- 
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FIG. 6: Shear-rate dependence of the off-diagonal element D^^^y of the thermal diffusion tensor for x\ = 0, K22 = ti2 and 
several values of the mass ratio /i = m\/m2. 



P 



5ik - 



(l-a*9a.) 



p* _ p* 



(4.18) 



mim2 \ mim2 / 

In order to get these expressions use has been made of the identity 

(fca + a)"^ = 6-i]l-6-2a, (4.20) 

where h is an arbitrary constant and a is the tensor with elements — aSixSjy. 

In the case of mechanically equivalent particles, Dp^ij — Dx^ij — and Eq. (|4.17p reduces to 



A, = I ku - ^7 , . (4.21) 

Equation (|4.2ip gives the self-diffusion tensor of tagged particles under thermostatted USFj^ For a general binary 
mixture, the expression (|4.17p for the diffusion tensor D^- coincides with the one derived beforei^^ in a stationary 
state with the constraints p = const and T = const. Finally, it is also apparent that, except for vanishing shear rates, 
the expressions of the generalized transport coefficients (|4.17p - (l4.19p in the thermostatted state differ from the ones 
derived in the absence of a thermostat, Eqs. (I4.7p - (l4.9p . This shows again that the presence of the thermostat affects 
the transport properties of the system. 



V. ILLUSTRATIVE EXAMPLES IN THE TRACER LIMIT 



The results obtained in the preceding Section give all the relevant information on the influence of shear flow on 
the mass transport. In general, the elements Dij^ ^p-ij and I?T,ij present a complex dependence on the shear rate 
and the parameters of the mixture without any restriction on their values. However, although the solution to Eqs. 
(I4.7p - (|4.9p (in the unthermostatted case) and Eqs. (|4.17p - (l4.19p (in the thermostatted case) is simple, it involves 

quite a tedious algebra due to the complex dependence of the partial pressure tensors P^^jj and the thermostat 
parameter a on the mole fraction xi and the reduced shear rate a* . To show the shear-rate dependence of the tensors 
Tij = {Dij, Dp^ij, Dx^ij} in a clearer way, the tracer limit (xi — > 0) will be considered here in detail. In addition, 
to make some contact with computer simulation results, the thermostatted case will be studied. In the tracer limit 
case, P ~ P2 and the partial pressure tensors Pi and P2 have a more simplified formSf^^ In particular, dxiPj;p — 

and -Pi = ^1%/^^- The exphcit expressions of the partial pressure tensors in the tracer limit are provided in 
Appendix [Cl 
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As expected, T^z = T^x = Tyz = T^y = 0, in agreement with the symmetry of the problem. As a consequence, 
there are five relevant elements: the three diagonal {T^x, Tyy, and T^z) and two off-diagonal elements {Txy and Tyx). 
In addition, Txx 7^ Tyy = Tzz and Txy ^ Tyx- The equality Ps.yy = Ps.zz implies Tyy = Tzz- This property is a 
consequence of the interaction model considered since for non-Maxwell molecules computer simulations show that the 
yy and zz elements of the pressure tensor are different. In Figs. [T]-[6l the relevant elements of tensors D*j, D*j^j 
and are plotted as functions of the reduced shear rate a* for K12 — K22 and several values of the mass ratio 11. 

Here, the tensors have been reduced with respect to their Navier-Stokes values (except Dx^ij), namely, D*j = Dij/Do, 
Dpij = Dp^ij/Dpfi and f^.y — Dr^ij/xiDo. One third of the trace of these tensors is plotted in Figs. [T]-[3l while 
the xy element is plotted in Figs. HHSl We observe that in general the influence of shear flow on the mass transport 
is quite important. It is also apparent that the anisotropy of the system, as measured by the traces ■^D'^f^, kk' 
and grows with the shear rate. This anisotropy is more significant when the impurity is heavier than the 

particles of the gas. Moreover, the shear field induces cross effects in the diffusion of particles. This is measured by 
the (reduced) off-diagonal elements D*y, D* ^y and D^^y. These coefficients give the transport of mass along the x 
axis due to gradients parallel to the y axis. While D^^ and D* ,^y are negative, the coefficient DT,xy can be positive 
in the region of small shear rates. We observe that, regardless of the mass ratio, the shapes of D*y and D* ^y are 
quite similar: there is a region of values of a* for which —D*y and —D* ^y increase with increasing shear rate, while 
the opposite happens for larger shear rates. The magnitude of T>^ is smaller than that of the elements —D*y and 
~D*^y, especially when the tracer particles are lighter than the particles of the gas. In this latter case, is 
practically negligible. 

VI. DISCUSSION 

Diffusion of particles in a binary mixture in non- Newtonian regimes is a subject of great interest from a fundamental 
and practical points of view. If the mixture is strongly sheared, the mass flux ji can be significantly affected by the 
presence of shear flow so that the corresponding transport coefficients may differ significantly from their equilibrium 
values. In addition, the resulting mass transport is anisotropic and thus it cannot be described by scalar transport 
coefficients but by shear-rate dependent tensorial quantities whose explicit determination has been the main objective 
of this paper. 

In order to gain some insight into this complex problem, a dilute binary mixture of Maxwell molecules under 
USF has been considered. This is perhaps the only interaction potential for which the Boltzmann equation can be 
exactly solved in some specific non-homogenous situations, such as in the case of the USF problem. In particular, the 
corresponding rheological properties of the mixture (nonlinear shear viscosity and viscometric functions) have been 
obtained for arbitrary values of the shear rate and without any restriction on the parameters of the mixture (masses, 
concentrations, and force constants).* This exact solution is of great significance in providing insight into the type 
of phenomena that can occur in conditions far away from equilibrium. In this paper, the interest has been focused 
on situations that slightly deviate from the USF by small spatial gradients. Under these conditions, a generalized 
Chapman-Enskog method^'^!-:^ around the shear flow distribution has been used to determine mass transport in the 
first order of the deviations of the hydrodynamic field gradients from their values in the reference shear flow state 
fs^\ In this case, the mass flux j^^"* is given by Eq. p.24p . where the corresponding set of generalized transport 
coefficients {Dij, Dp^ij, DT.ij} are the solutions of the coupled algebraic equations (|4.7p - (|4.9p in the unthermostatted 
case, while they are explicitly given by Eqs. (|4.17p - (|4.19[) in the presence of a Gaussian thermostat. This type of 
external forces are usually employed in nonequilibrium molecular dynamics simulations to compensate exactly for the 
viscous increase of temperature. 

As expected, the results show that the coefficients {Dij, Dp^ij , Dx^ij} present a complex dependence on the shear 
rate and on the masses, mole fractions, and force constants. This is clearly illustrated in Figs. [THB] for the tracer 
limit case (xi — *■ 0). The deviations of {Dij,Dp^ij,DT^ij} from their equilibrium values are basically due to three 
different reasons. First, the presence of shear flow modifies the collision frequency of the conventional diffusion problem 
(/oAi2/miTO2) by a shear-rate dependent term. Second, given that the binary mixture is in general constituted by 

particles mechanically different, the reference shear flow states z}"' and /j"'' are completely different. This effect gives 
rise to terms proportional to Pj^^^ — {pi/p)P*j. Third, in the unthermostatted case, the generahzed coefficients are 
coupled due to the inherent non-Newtonian features of the USF state. Each one of the three effects is a different 
refiection of the extreme nonequilibrium conditions present in the mixture. 

It is apparent that the results presented here in the particular case of Maxwell molecules may be relevant for 
interpreting computer simulation results. Sarman, Evans, and Baranyai^* carried out time ago molecular dynamics 
simulations in a strongly sheared Lennard- Jones binary mixture to evaluate the self- and mutual-diffusion tensor by 
means of Green-Kubo formulae. To the best of my knowledge, this is the only computer experiment in which the 



13 



shear-rate dependence of the diffusion tensor Dij has been measured. They considered an equimolar Lennard- Jones 
mixture at two different densities and the parameters in the potential were adjusted to model an argon-krypton 
mixture, which means that the two components are fairly similar. As already said in Ref. [isl . when one considers 
this type of mixture (xi = 0.5, mi/m2 ~ 0.48, kh = K22 — K12) in the thermostatted case, the general qualitative 
dependence of the (reduced) mutual diffusion tensor Dij(a*)/ Dq on the (reduced) shear rate agrees quite well with 
computer simulations. Thus, theory and simulation predict that in general, the xx element increases to a maximum 
and then it decreases again, while the yy element decreases with increasing shear rate. The off-diagonal elements xy 
and yx are negative and their magnitude increases with a* for not very large values of the shear rate. However, kinetic 
theory predicts that \Dxy\ > \Dyx\, while the opposite happens in computer simulations. On the other hand, at a 
quantitative level, the influence of shear flow on diffusion is much more modest in the molecular dynamics simulations 
than the one found theoretically for dilute gases. This is probably due to the fact that the shear rates (in reduced 
units) applied in the simulations are not large enough to observe significant changes of the diffusion tensor relative 
to its equilibrium value. An alternative to overcome the difficulties for reaching large shear rates in nonequilibrium 
molecular dynamics at low-density is the direct simulation Monte Carlo method.— I hope that the results derived here 
for Maxwell molecules for Dij, Dp,ij and DT,ij stimulate the performance of Monte Carlo simulations to assess the 
reliability of the Maxwell results to describe mass transport in strongly sheared mixtures for more realistic interaction 
potentials. 

As said before, it must noted that, in order to observe large effects of shear flow on the tensors {Dij, Dp^ij, DT,ij}, 
the (reduced) shear rate must be at least of the order of 1. This means that for the inert gas fluids considered in 
this paper, non-Newtonian effects on mass transport could be observable for shear rates practically unattainable in 
the laboratory. In this sense, one should look at fluids that are observed to be non-Newtonian, such as colloidal 
suspensions, polymeric hquids, gels,- • • ^L2^ 

Although the results derived in this paper have been focused on the mass transport, the remaining transport 
coefficients associated with the pressure tensor P^p and the heat flux vector q'^"'^^ could be determined from the 
integral equations (|3.18p - (|3.2ip . Nevertheless, in practice this calculation cannot be carried out analytically by using 
the Boltzmann equation since the fourth-degree moments of USF (whose explicit expressions are not known in the 
Boltzmann equation, except for a single gas^) are needed to get the heat flux. In order to overcome such a difficulty 
one can use a model kinetic equation that preserves the essential features of the true Boltzmann equation but admits 
a more practical analysis. Perhaps the most well-known model for gas mixtures is the Gross-Krook (GK) kinetic 
model^S. In this model the Boltzmann operator Jrs[fr, fs] is replaced by the relaxation term 



Jrs[fr,fs] -Vrsifr " frs), (6.1) 



where 



and 



3/2 



2kBTr 



-(v - u^^)" 



(6.2) 



rrirVLr + m,u, 

Irs = , 6.3 

rrir + TTls 



rrirms 



{Ts~Tr) + -^{Mr-Visf 



(6.4) 



The partial temperatures are defined by Eq. 
given by 



For Maxwell molecules, the effective collision frequency i^rs is 



rrir + 
rrirms 



1/2 



(6.5) 



where A is a constant to be fixed by requiring that the model reproduces some transport coefficient of the Boltzmann 
equation. An exact solution to the GK kinetic model for a binary mixture in USF has been found. '^-'^ The comparison 
of the GK results with those from the Boltzmann equation* at the level of the rhcological properties shows good 
agreement, confirming the reliability of the GK model in computing transport properties in far from equilibrium 
situations as well. Starting from the USF solution of the GK model, the fiuxes P^^j^ and q*^^-' are obtained in 
Appendix iP] in the thermostatted case. With all the transport coefficients known, the constitutive equations for the 
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mass, momentum and heat fluxes are completed and the corresponding set of closed hydrodynamic equations for the 
mixture can be derived. This allows one to perform a linear stability analysis of the hydrodynamic equations with 
respect to the USF and determine the conditions for instabilities at long wavelengths. Previous results for a single 
gasS. have shown that USF is unstable when the perturbations are along the velocity gradient (y direction). The 
problem now is to extend this analysis to the case of multicomponent systems. Work along this line will be reported 
in the near future. 
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APPENDIX A: CHAPMAN-ENSKOG-LIKE EXPANSION 

In this Appendix, some technical details on the determination of the first-order approximation f[^^ by means of the 
Chapman-Enskog-like expansion are provided. Inserting the expansions (|3.9p and (|3.11[) into Eq. (|3.1a[) . one gets the 
kinetic equation for /}^^ 

di"'A'' - ^ {a.,V, + aV.) A« +a6u.^+ C,fi'^ + M,fi'^ = - [a^^ + (V + Uo) • v] f['\ (Al) 

The velocity dependence on the right-hand side of Eq. (|A1[) can be obtained from the macroscopic balance equations 
p.2p - p.4p to first order in the gradients. They arc given by 

di^^xi = -(uo + Su) ■ Vxi, (A2) 
di^^Su = -(uo + Su) ■ VSu - p-^V ■ P(°), (A3) 

dl^^p = -(uo + du) -Vp-pV -Su - - i^aPg^ + P(") : V6uj , (A4) 

di^^T - -(uo + 5u) ■ VT - |- (aPiJ) + p(o) : VJu) , (A5) 
where use has been made of the result j^"-* = q'"' = 0. In addition, 

Use of Eqs. in Eq. ^M^ yields 

+Bi -Vp + Ci-yr +Di: VSu, (A7) 



where 



AiAc) ^ — T~ci + -—^ A8) 

oxi p oouj oxi 



i3M(c) = -^c, + i^^, (A9) 
op p oouj op 
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Upon writing Eq. (jAlip use has been made of the expression of the total pressure tensor P^j ' of the mixturo^i^ 

P<« . (A12) 

where ?7yfef is the viscosity tensor. 

The solution to Eq. (jA7p has the form given by Eq. (|3.17p . where the coefficients A-i, Bi, Ci, and Vi arc functions 
of the peculiar velocity and the hydrodynamic fields xi, p, T, and Su. The time derivative acting on these quantities 
can be evaluated with the replacement 

di°^ - - + 2«) ipdp + Tdr) . (A13) 

Moreover, there are contributions from df*^ acting on the pressure, temperature, and velocity gradients given by 

I 3 dxi dxi j dp dp j 

-h-'-^^2p'AvT, (A14) 

9fVT = -v(^aPi°}+2Ta 



3p 



xy 



\ 3p dxi dxi I \ 3p2 3p dp dp I 

df^^ViSuj = Vidf^Suj = —Gjk'ViSuk- (A16) 

The corresponding integral equations (|3.18|l - (|3.20p can be obtained when one identifies coefficients of independent 
gradients in Eq. (|A7[) and takes into account Eqs. (jA14[) - (|A16p and the mathematical property 

|aP<;.+2a)(p| + r^)x + a„*„,f. (A17) 

where in the last step it has been taken into account that X depends on Su through c = V — Su. 

APPENDIX B: GENERALIZED TRANSPORT COEFFICIENTS ASSOCIATED WITH THE MASS 

TRANSPORT 

In the unthermostatted case (a = 0), the integral equations defining the generalized transport coefficients Dij, Dp ij 
and DT,ij are given by Eqs. (|4.ip - (|4.3p . To get these coefficients, one multiplies (|4.iP " (|4.3p by miCi and integrates 
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over velocity. The result is 



aP(,lHpdp + TdT) 



3p 



mi?7i2n 



P 



mim2n ( pXi2 



mim2 



nil / dc CiAi, 



2ap 
~3f 



-aP!^^Jipdp + TdT){^D. 



aikDp,kj + 



Phil. _ ^^d^pio}] u 



mi 7712 3p 



7771 J dc C,B,,j + ^DT.^Ja - Pdp)P^^l\ 

aP^'J ipdp + Tdr) {^Dr,,) - f [a.kDr^k, + - |(1 + Tdr)P^oj^ Dr,, 

= 7771 Idc c.Ci,, - ^Dp^^.idrPil^), 



2 

3p 



where 



J dc c,c,/i°). 

Upon writing Eqs. (|Bl[) - (jB3[) . use has been made of the relation (|2.12p . which yields the results 

7771 J dc a (£1^1 + A^1^2) = ~n\l2Dij, 

dc Q (£iBi + M1B2) = —Dp ij, 

777 1 7772P 



7771 J dc Ci {C1C1+M1C2 



P^Ai2 



D 



T,ij ■ 



77717712^ 

The velocity integrals appearing in Eqs. (|Bip - (jB3p can be performed by using Eqs. 



a ( p(0) _ Pi p(o) 

'"P \ ^1,13 p ^i] 



7^1 j dc c,Ai,, = - [d,M% - jd^^P^P) ' 
nil J dc CiJ3ij 
mi J dc CiCij 



f, I p(0) Pi p(0) 



(Bl) 



(B2) 



(B3) 
(B4) 

(B5a) 
(B5b) 
(B5c) 

(B6) 
(B7) 
(B8) 



The generalized transport coefficients Dij, Dp^ij, and Dx^ij can be written as Dij — Z)o-D*j(a*), Dp^ij — Dp ^D* j^ (a*), 
and Dr.ij — Dt.qD'^ ^^{a*) where , Dp^^, and D^^^ are dimensionless functions of the shear rate. Moreover, from 
dimensional analysis, Da ^ T, Dp^ ^ T'^ /p, and Dt.o ^ /p. Therefore, 



ipdp + Tdr) ( ^^D^, ] = ipdp + Tdr) { ^^DoDf^ ] = ™^D, 



{pdp + Tdr) { '-Dp,^, 



{pdp + Tdr) [ ^Dp,oD;^^^ ) = 0, 



(pdp + Tdr) {^Dr,^,) = {pdp+Tdr) [^Dr,oD*r,,,) = 0, 
where use has been made of the identity 

{pdp + Tdr) X{a*) = [da'X) {pdpU* + Tdra*) - 0, 
with a* — a/Q ^ T/p. Taking into account the above results one arrives at the set of algebraic equations (I4.7p - (|4.9p . 



(B9) 
(BIO) 
(Bll) 
(B12) 
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APPENDIX C: RHEOLOGICAL PROPERTIES IN THE TRACER LIMIT 

The explicit expressions for the pressure tensors P*^j = Pg^ij/p in the USF are provided in this Appendix for the 
special case of tracer limit (a;i — > 0). The nonzero elements of P2 ij are given hy^ 

Plyy = Plzz = 1 + 2L0a* ' ^^^^ 

_ 1 + 6u;a* 

- 1 + 2u;a* ' ' 
P* - -3— - (c-i) 



where a* = a/C, a* — a/C,, ( being defined by Eq. (|4.4p . Moreover, 



2 /K22I + M /^.x 
722 = W , (C4) 



722(1 + A*) ' V K12 M 

where /i = mi/m2 is the mass ratio. The (reduced) thermostat parameter is given by a* — max(Q!o, ap) where^i^ 

"0 = -^G{uja*), ofg = ^G{2p,a*) - ^712, (C5) 

where G{z) — | sinh^[i cosh^^(l + 9z^)] and 712 — \i2/^'i2 — 0.648. Usually, uq > except for very large shear 
rates and/or very disparate mass binary mixtures 1^ 



The nonzero elements of Pi ij are given hy^ 



+ i(l + 2a;a*)(2a*+e)2|, (C6) 

" A(l + La*) {^^'' ~ + + ^^^^^^ ~ + + ^"""^^ (^"* + ' ~ 

+ i(l + 2Loa*) [{2a* + ef + 2a*^) + (2712 - l)a*'] | , (C7) 

r>* xicja* f, * * \ / l,l + 2ci;a* 

PU = -TT^^T-^ (712 - 7T)(2a* + 7i2)(2a* + e) + (712 - 7.)" 



A(l + 2tja*)2 [^'"-^ 2'^ ' 2' u 

X (2a* + e + 2uoa*f3) + —(1 + 2uja* f (2a* + e) I , (C8) 
2tj 



where e = 712 = 1/2/i, and 

A=(2a* + e)2(2a*+7i2)-^/3a*2. (C9) 
APPENDIX D: MOMENTUM AND HEAT TRANSPORT AROUND USF FROM THE GK MODEL 



model (|6.ip . The first order corrections to the fluxes are 



This Appendix addresses the evaluation of the fluxes P^^ and q*^^^ in the thermosttated case from the GK kinetic 
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where the partial contributions to the transport coefhcients are defined as 



-nis j dc CtCjVs^keic), (D3) 
Jdcc'c^AsAc), (D4) 
,y = - ^ / dc c2c,6,j (c), (D5) 



J acc^c^C.^c). (D6) 



From the above partial contributions one can get the generalized shear viscosity rjijke = 'r]i.ijke + fl2.ijke, the generalized 
Duffour coefficient D'/j = ij + ^'iij^ the generalized pressure energy coefficient = Li^ij + L2.ij^ and the 
generalized thermal conductivity Ay = Ai^y + \2,ij- The quantities {yls^i, Cs,i, X's^ij} still verify the integral 
equations p.l8p - p.2ip (with a — —aP^/Zp) with the only replacement 



Ci.f['^ + Mif^'^ ^ v,,f['^ - - i.i2/l^\ (D7) 



where vi = vn + 1^12, 



and 



AO) AO) 
f(l) _ /ll „ :(1) _ ■/12 / „ V ,(1) /p>o^ 

fio) = ( y/^^p J!VC^) . (D9) 

In Eq. (IDSp . Urs — mr/{mr + TOs)- 

In order to get the coefficients rjs^ijke, D'^iji Ls,ij and Xs.ij, it is convenient to introduce the velocity moments 

4:1,™ = / dcc^^CA., (DlO) 
nx™ = / dcc^4CSi,„ (Dll) 
43,™ = / dc c^^CCi,., (D12) 

<ll = / dc c^^cr^M,, (D13) 

and similar definitions for the species 2. The knowledge of the above moments allows one to get the expressions of 
the coefficients rii,ijke, , and Ai,.y . The method to evaluate the moments xj^^^^^^, Y^^ll^, -^^^l ^, and W^^^^ 

(i) 

is quite similar. Here, as an example, the mathematical steps to determine the moments Xj^ ^ ^ associated with the 
transport coefficients D'{ will be analyzed in detail. First, in the thermosttated case, Eq. (|3.18p with the change 
(|D7P becomes 

f ^ , ^ \ A , TOi?7i2n / i^ii (0) 1^12 I \A0)\ n /I ^T^^A\ 

- ^Cy- 1/1+ a— ■ c Al,^ + , ^ /ll + 7r^Aii2 ?^2 - "1 /12 CjDji = Ai^„ (D14) 

V dc I p VnifesTi nin2fc_BTi2 / 
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where Ai is given by Eq. (jASp . Upon writing (|D14p use has been made of the constitutive form p.24p for the mass 
flux. Now, we multiply Eq. (|D14p by c^c^c™ and integrate over velocity. After some algebra, we get 



a^4^2i^,+i^„, + Wi + ik + e + m)a] Xl[ 



(0 



R 



ii) 

k,£,m ' 



(D15) 



where 



-J- TOir7i2n / 2kBT 

■^i.i 

P \ mi 



(k+e+m+l)/2 



Vll (fe+£+m-l)/2 



Vl2 , N (fc+f+T?x-l)/2 



n2kBT 



(D16) 



/d 



+ - {SjxkMk-l,£,7n + SjylMkJ-l.m + 5j^mMk,l,m-l) ■ (D17) 

In Eqs. (|D16p and (jD17p . we have introduced the temperature ratios xi = Ti/T and Xi2 = T12/T and the (unper- 
turbed) moments of the USE 



= J dcc':clc^fi°\c). (D18) 



The explicit shear-rate dependence of xi, X12 and Mk.e.m can be found in Ref. l3ll Moreover, 

AM,™---^/^r(^)r(i±l)r(!^) (di9) 

if {k,£,m) are even, being zero otherwise. The solution to Eq. (|D15p can be written as 

^M,™ = E (fc^(-«)'[^i + {k + i + . (D20) 

Note that Eq. (|D20p is still formal since one needs to know the coefficients Dij . They can be consistently determined 
from their definitions (|3.25p . Once these coefficients are known, Eq. (|D20p allows one to get the coefficients D'^ ij- 

The same method can be applied to evaluate the remaining moments. The moments Y^^j and Z^*^ are given by 

^S™ = E + (fc+^+-)"r^^+'^5«,,,+,,^ , (D21) 



where 



fc! 



{~ay[iyi + ik + £ + m)a] T^!^ 



(D22) 



^k,e,m 



— Bi,i — 



P f^ksT 



(fe+£+m+l)/2 



p \ mi 

X Afc+iSjx ■e+Sjy ,rn+Sj^ Dp,3i ) 



1^11 (k+l+Tn-l)/2 , 1^12 , N (fe+f-|-m-l)/2 

T^Xi H ; — ;7;Mi2l"2 - '^ijXi2 

ksT n2kBT 



(D23) 



(0 



k.l.m 



(/c+£+m+l)/2 r 



J \ mi 

X /+<5js ,m+(5j^ DT,ji ■ 



Vll (k+l+m-l)/2 , V12 

ksT n2kBT 



I N (fe+f+m-l)/2 



(D24) 
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The expressions of Bi^i and Ci^i are formally identical to that of Ai i, except that the operator dx^ appearing in 
(|D17p must be replaced by the operators dp and dx in the cases of Bi i and Ci^i, respectively. Finally, the expression 



where 



k,i,m 



(y) 



fc! 



^V" k,e,m. 



d 



(0) 



^ dp 



d 

T — 

dT 



Mk-1 ,£,m+l } 
-mSiz {SjxMk+l/'m~l + 5jyMkA+l,rn-l) ■ 



(D25) 



(D26) 



The transport coefficients D'^^j, Li ^j, Xi^ij, and rji^ijki can be obtained from Eqs. (|D20p . (jD2ip . (|D22p and (|D25p . 
respectively, in terms of the shear rate and the parameters of the mixture. Their respective counterparts for species 
2 can be easily determined from them by making the changes: toi m2, ni n2, and kh K22- The expression 
of the Duffour tensor D'^'j coincides with the one obtained before^'* in a stationary state with Vp — VT — 0. Finally, 
note that the USF moments Mh p „ can be written as^ 



ni 



2kBT 
nil 



(k+i+m)/2 



M, 



where the dimensionless moments M, 



k,i,; 



pdpMk/,m = ni 



depend on p and T through their dependence on a* . Consequently, 
2kBT'^ (fc+''+™)/2 



mi 



il-a*da.)M*,, 



(D27) 
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+ a* da. 
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